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ABSTRACT 



Stiffened plates are designed for minimum volume si±)ject 
to constraints on: Von Mises maximum stresses, nodal dis- 

placements, height to thickness ratios of frame webs, and 
width to thickness ratios of frame flanges. Design variables 
are plate thicknesses and stiffener dimensions. 

A finite element analysis program is developed for the 
design of stiffened plates using numerical optimization 
techniques. The program may be used as a stand alone analy- 
sis tool or may be coupled to an optimizer of user's choice. 
Rectangular plate elements and frame elements are used for 
the idealization of stiffened plates. 

Design examples are presented tC' demonstrate the design 
method. 
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INTRODUCTION 



The development of high-! 5 peed digital computers has made 
possible significant changes in the structure design process. 
One of these is the availability of various mathematical 
programming methods for use in design optimization. The 
computer's speed allows the designer to now consider a much 
wider range of design alternatives. The optimization proce~ 
dure provides a means of systematically choosing from among 
these alternatives based on some predetermined rational 
criterion. 

Even when the selected numerical method is able to arrive 
at the optimum design, tne result is only as good as the 
design model. Here there is an even greater need for design 
experience and sound engineering judgment. The design model 
must be carefully developed to realistically represent the 
design in question. 

The finite element model for stiffened plate is shown in 
Fig. 1.1, where the plate has been descritized by rectangular 
plate elements and the stiffener by frame type elements. The 
eccentricity of stiffener is transformed to the linked nodal 
point by applying a linear equation that realtes displacement 
degrees of freedom. 

The purpose of this thesis is to develop a finite element 
analysis program for stiffened plates, and to design the 
optimum stiffened plate by coupling two programs; the 
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Figure 1.1 Stiffened Plate Model 
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analysis program and a general purpose non-linear optimizer, 
COPES/CONMIN [Refs. 1-2]. 

The remaining chapters of this thesis are outlined as 
follows ; 

Chapter II presents briefly the finite element method 
used in the analysis program for stiffened plates. 

Chapter III presents the basic concepts of the cptiiri- 
zation methods used in the COPES/CONMIN. 

Chapter IV presents design examples . 

Chapter V offers conclusions and recommendations. 
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II . 



FINITE ELEMENT METHOD 



A. GENERAL 

The finite element method is now well established as an 
engineering tool of wide application. The fundamental 
concept of the finite element raethod is that any continuous 
quantity, such as pressure or displacements, can be approxi- 
mated by a discrete model composed of a set of piecewise 
continuous functions defined over a finite number of si±»- 
domains. The piecewise contin\ious fxanctions are defined 
using the values of the continuous quantity at a finite 
number of points in its domain,. 

The formulation of the fin..te element method can be traced 
to energy procedures, principally the minimum potential energy 
method and the minimiam complimt2ntary energy method. The 
minimum potential energy method is associated with assumed 
displacement parameters as unknowns and is usually termed 
the "displacement" or "stif fnetss" method. On the other hand, 
the minimum complimentary energy method dealt with stress 
parameters and is termed the "flexibility" or "force" approach. 
The ease with which a continuous displacement pattern can be 
prescribed (compared to the alternative approach of forming 
an equilibrating internal force field) has aided the wide- 
spread use and development of the finite element displace- 
ment approach. The displacement model and the stiffness 
analysis are employed in the analysis program developed here. 
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This chapter will briefly present some of the general 
concepts of the finite element method used in the analysis 
program. 

B. FINITE ELEMENT DISPLACEMENT APPROACH 

The displacement formulation involves derivation of the 
stiffness matrix of each individual element. The stiffness 
matrix of the entire assembled structure is then obtained by 
super-position. This matrix, along with the prescribed 
displacement boundary conditions and loads, is used for the 
solution of displacements and stresses. 

1 . Element Analysis 

For the structural applications at least, the 
governing equilibrium equations can be obtained by mini- 
mizing the total potential energy of the system. The total 
potential energy, tt can be expressed as 

tr = i / £ dV - J u*^ f dV - J u*^ q dS (2.1) 

Zv-- v~'' S~~ 

where a and £ are the stress and strain vectors respectively, 
u the displacements at any point, f the body forces per unit 
volume and q the applied surface tractions. Integrations are 
taken over the volume V of the structure and loaded surface 
area S. 

The first term on the right hand side of Eq. (2.1) 
represents the internal strain energy and the second and 
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third terms are respectively the work contributions of the 
body forces and distributed surface loads. 

In the finite element displacement method, the basic 
steps for derivation of the element stiffness matrix are: 

a. Express the displacements to have unknown values only 
at the nodal points, so that the variation within any 
element is described in terms of the nodal values by 
means of interpolation functions. Thus 

u = N u^ (2.2) 

where N is the set of interpolation functions termed 
the shape functions and u^ is the vector of nodal 
displacements of the element. 

b. Express the strains within the element from the 
element nodal displacements as 

e = B u® (2.3) 

where B is the strain-displacement matrix generally 
composed of derivatives of the shape functions. 

c. Express the stresses relating to the strains by use of 
an elasticity matrix D containing the appropriate 
material properties, as follows: 

a = D £ (2.4) 

d. Establish the equilibrium equation of element. 
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Provided that the element shape functions have been chosen 
so that no singularities exist in the integrands of the 
functional, the total potential energy of the continuum will 
be the sum of the energy contributions of the individual 
elements. Thus 



TT = J 1T^ (2.5) 

e 

where ir^ represents the total potential energy of element e 
which, on use of Eq. (2.1), can be written 



= j f D B u® dV - / u®"^ 



f dV 



( eT T 
- J N q dS 



( 2 . 6 ) 



where V is the element volume and S the loaded element 
e e 

surface area. Performance of the minimization for the 
element e with respect to the nodal displacements u® f or the 
element results in 




/ (B*^DB) u^ dV - / f dV - / q dS 

V v = S'"'' 

e e 



= k'= 



u 



- P 



(2.7) 



where 
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( 2 . 8 ) 



e 

P 




f dV + / 

s ^ 

e 



q dS 
1 . e 



are the equivalent nodal forces for the element, and 



k® = / D B dV (2.9) 

= V = = = ® 

e 



is termed the element stiffness matrix. 

2 . Direct Stiffness Method 

The real elastic structure is now represented by a 
finite n\amber of small, discrete elements. Once their 
approximate behaviors, identified by their individual stiff- 
ness matrices k®of Eq. (2.9), have been established, the 
stiffness matrix K for the complete structure is obtained by 
the proper summation of each element stiffness matrix in the 
structure. The suinmation of the terms in Eq. (2.7) over all 
the elements, when equated to zero, results in a system of 
equilibrium equations for the complete continuum. This 
assembly process is known as the Direct Stiffness Method. 
These equations are then solved by any standard technique to 
yield the nodal displacements. Note that K is symmetric and 
positive-definite. 

C. FINITE ELEMENTS USED 

The finite elements used in the analysis program will be 
described briefly in this section. 
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1. Rectangular Plate Element 



The rectangular plate element used here is illus- 
trated in Fig. 2.1, where each nodal point has 6 degrees of 
freedom, 3 transitional displacements and 3 rotational 
displacements. This element with 4 corner nodal points has 
the element stiffness matrix of order 24 x 24. The corres- 
ponding displacements of each node u^^ , U 2 / . . . , U 2 ^ will be 
taken to be positive in positive directions of the xyz- 
coordinates . 

The main assumptions in the method are that displace- 
ments are small compared with plate thickness, the stress 
normal to the midsurface of the plate is negligible, and 
normals to the midsurface before deformation remain straight 
but not necessarily normal to the midsurface after deformation. 

The assumed displacement functions will be taken to 

be linearly varying in the plane of the element. These 

displacement functions will ensure both deflection and slope 

compatibility of the adjacent elements. The stiffness matrix 

of the rectangular plate element, which is matrix product 
T 

B DB of Eq. (2.9) integrated over its volume, is summarized 
by Przemieniecki [Ref. 3]. 

Note that element stiffness matrices are formed 
directly in the global coordinate system so that no trans- 
formations from local to global coordinates are required. 

2 . Frame Element 

The frame element as a stiffener has 6 degrees of 
freedom for each node, such as those of plate elements. The 
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2, z' LOCAL AND GLOBAL 
COORDINATE AXES 



Figure 2.1 Rectangular Plate Element 
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basic assumptions are that the stiffener element is straight 
with uniform cross section capable of resisting axial forces, 
bending moments about the two principal axes in the plane of 
its cross-section, and twisting moments about centroidal axis, 
and that its deflection due to shearing strains are neglected. 
The stiffener element has 2 kinds of options which are illus- 
trated in Fig. 2.3. The width and thickness of flange of 
the 3rectangular stiffener element may be referred as zeros. 

In order to determine the stiffness property of a 
complete structiore, a global coordinate must be established 
for all xinassembled structural elements so that all the 
displacements and their corresponding forces will be referred 
to this system. Since the element stiffness matrices k are 
initially calculated in local coordinates, suitably oriented 
to minimize the computing effort, it is necessary to intro- 
duce transformation matrices changing the frame of reference 
from a local to a global coordinate system. The first step 
in deriving such a transformation is to obtain a matrix rela- 
tionship between the element displacements u in the local 
system and the element displacements u* in the global system. 
This relationship is expressed by the matrix equation 

u = T u' (2.10) 

where T is a matrix of coefficients obtained by resolving 
global displacements in the directions of local coordinates. 

The transformation matrix T is given by 
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Figure 2.2 Frame Element 
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a. Rectangular type frame stiffener 




b. T type frame stiffener 




Figure 2.3 Stiffener Options 
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represents matrices of direction cosines for local x, y and 
2 directions, respectively, measured in global system x, y 
and z. Similarly the relationship for forces is expressed 
by the: matrix equation 



p = T p ' (2.11) 

where p is the force vector in the local system and p' is 

the force vector in the global system. Matrix A is orthogonal; 

-1 T 

that is, T = T . Therefore 

u' = u (2.12) 

and 

p' = t'^p (2.13) 
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Let u' and u be two ways to describe the same virtual dis- 
placement. Virtual work is 




6u p 



T T 

6u' T p 



(2.14) 



So 



6u'‘^(p' - t'^P) = 0, p' = T*^p (2.15) 

To transform the stiffness matrix, we start with k'u' = p' 
and substitute from the preceding equations. 

k'u' = p' = t"^p = T"^k u = T"^k Tu ' (2.16) 

Since this relation is presumed valid for any u', we conclude 
that the required stiffness transformation is 

k' = k T (2.17) 



D. OFFSETTING OF RIGID LINKS 

One of the most important advantages of the finite ele- 
ment technique is that an assembly of different structure 
elements such as plates and frames can be dealt within a 
single coordinate analysis. Usually the neutral surfaces of 
plate and stiffener are not coincident: the stiffener is on 

one side of the plate. A standard preliminary treatment is 
to connect adjacent plate and stiffener nodes by a rigid 
link, so that degrees of freedom of the stiffener are replaced 
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by degrees of freedom of the plate. The usual assembly is 
then possible. The necessary transformation is now described. 

The stiffener element in Fig. 2.4 has the usual 12 degrees 
of freedom--6 at node A and 6 at node E. With reference to 
these degrees of freedom, element load and stiffness matrices 
are p and k. Similar degrees of freedom are used at nodes 1 
and 2 of plate element of rigid links A-1 and B-2. The 
"master" degrees of freedom at node 1 and "slave" degrees 
of freedom at node A have the relation 
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Figure 2.4 Rigid Links 
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A similar expression is written for link B-2 by replacing 
subscripts A and 1 by B and 2. The transformation vectors 
p and matrix k, associated with degrees of freedom at nodes 1 
and 2 , are 



p = T p 



(2.19) 



T — 

k = T k T 



( 2 . 20 ) 



where 



?12xl2 



A, 0 






The foregoing transformation makes the translational 
displacements depend on the rotational displacements and so 
introduces an unwanted quadratic field into the translational 
displacements. It is found that a frame-stiffened plate is 
overly flexible in a coarse mesh. Mesh refinement helps: 
error is reduced by a factor of four if the number of elements 
is doubled. The error can be eliminated through the addition 
of one more degree of freedom. Further details of the error 
reduction are described in [Ref. 4], 
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III. OPTIMIZATION TECHNIQUE S 



A. INTRODUCTION 

In this chapter some of the fundamental ideas and formu- 
lation methods of mathematical programming are introduced 
to understand the optimization method and the optimizer 
COPES/CONMIN. Design parameters used here are then dis- 
cussed. Fox [Ref. 5] and Himmelblau [Ref. 6] provide an 
extensive discussion of numerical optimization techniques 
and their application to engineering design. 

B. DEFINITIONS 

In discussing the optimization methods, the following 
definitions will be useful: 

Design variables — the design variables are the numerical 
parameters for which values are to be chosen in producing a 
design. In a structural problem, they might be plate thick- 
ness, frame dimensions, etc. 

Objective function — the objective function, is the single 
valued function with respect to which the design is optimized. 
In a structural design problem, it might be the weight, volume 
or fabricated cost of the structure. The selection of an 
objective function can be one of the most important decisions 
in the optimum design process. 

Constraints — practical design problems are usually sub- 
ject to a series of constraints which must be satisfied in 
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order to produce an acceptable design. These constraints 
may be linear or non-linear. If a parameter is beyond the 
value of a specified value, it is said to be violated. 

Side constraint — the side constraint is a constraint 
which restricts the specified range of a design variable 
for reasons other than the direct consideration of performance. 

Feasible design — the feasible design (or acceptable 
design) is a design which satisfies the specified constraints. 

Infeasible design — the infeasible design (or unaacep- 
table design) is a design which does not satisfy the speci- 
fied constraints. 

The general nonlinear optimization problem is expressed 
mathematically as follows: 



Minimize F(x) (3.1) 

Subject to: 

g^ (x) £ 0 for j = l,m (3.2) 

x^ < X. < xV for i = l,n (3.3) 

1 — 1—1 



where the vector, x, is the vector of n design variables. 

The objective function, F (x) , given by Eq. (3.1), as well 
as the constraint functions given by Eq. (3.2), may be linear 
or non-linear functions of the design variables. They may be 
explicit or implicit functions of x, but must have continuous 
first derivatives. If it is desired to maximize F(x), the 
minimization of the objective function is used since maximum 
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F(x) can be treated minimum of -F (x) . The function (x) 
is the set of inequality constraints to be met. The m in- 
equality constraints must be satisfied to be a feasible 

LI 

design. Side constraxnts, x^ and x^, are the lower and upper 
limitations placed on the design variables. Side constraints 
could be included in Eq. (3.2), but are treated separately 
for efficiency. 

The n-dimensional space spanned by the design variables 
x^ is referred to as the design space. As stated previously, 
any design which satisfies the inequalities of Eq. (3.2) is 
referred to as a feasible design. If the design violates one 
more of the inequalities, it is said to be infeasible. The 
minimum feasible design is said to be optimal. 

Most nonlinear optimization programs update the vector 
of design variables by the iterative relationship 



X 



q+1 _ 



= + a= 



(3.4) 



where q is the iterative number, vector S is the direction 
of search in the design space, and a* is a scalar which 
defines the distance of travel in the direction S during the 
qth iteration. An initial design defined by X° must be 
defined and may be a feasible or infeasible design. 

The optimization process then proceeds in two steps. The 
first is the finding of S which will improve the design 
without violating constraints and the second is the determina- 
tion of a* which will improve the design as much as possible 
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when moving in this direction. The process is then repeated 
until there is no further design improvement, indicating 
that this is the optimum attainable design. 

C . COPES/CONMIN 

The COPES/CONMIN optimization program is a general pur- 
pose, non-linear optimizer capable of handling large, con- 
strained problems. This includes the conjugate direction 
method of Fletcher and Reeves [Ref. 7] for unconstrained 
function minimization and a modification of Zoutendijk's 
Method [Ref. 8] of Feasible Directions for constrained func- 
tion minimization. It has been successfully used in connec- 
tion with structural optimization [Ref. 9], airfoil design 
[Ref. 10], aircraft synthesis [Ref. 11], ani numerous other 
engineering applications. 

It was necessary to develop a subroutine, ANALIZ, which 
for a given design, would analyze a stiffened plate, and which 
would be suitable for coupling with the optimizer. The 
common block GLOBCM is required to couple the analysis sub- 
routine directly to the COPES/CONMIN. All variables, which 
are used as objective function, constraints and design varia- 
bles, must be listed in the common statement and the statement 
must appear in each subroutine the variables are used in. 

It is used by the optimizer as a catalog to identify where 
the design variables, objective function and constraints are, 
and what purpose they fulfill. 
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In order to execute the COPES/CONMIN program it is 
necessary to provide formatted data for COPES/CONMIN, followed 
by data for the ANALIZ program. Further details for the 
COPES data are explained in [Ref. 1] . 

D. DESIGN PAFAMETERS 

The design variables, constraints and objective function 
used in the dessign process are discussed in this section. 

1 . Design Variables 

The design variables are plate thicknesses, t^, and 
stiffener dimtmsions, being the web heights h^, the web 
thicknesses the flange widths W^, and the flange thick- 
nesses i = l,ns, where ns is the number of the different 

stiffener dimcjnsion sets. Thus, the total number of design 
variables depends on the sets of plate thicknesses and the 
sets of different stiffener dimensions. 

2 . Objective Function 

Total structure volume is considered as the objec- 
tive function in the design process. 

m ns 

Voliime = y V . + y N. V . 

i=l j£l ^ 

where is the volume of the ith plate element, m is the 

nimber of plate elements, is the volume of the jth stiffener 

set, and Nj is the n;mber of stiffeners in this. The number 
of stiffener sets is ns. 
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3. Constraints 



Design constraints are Von Mises maximum stress, 
nodal displacements, height to thickness ratios of frame 
webs, and width to thickness ratios of frame flanges. 
Stress: 



a 



a 



^ / 

. ./a 

13 



^ / 



max 



max 



1 < 0 



1 < 0 



i = l,ne 

j = 



i = l,ne 

j = 



u 2 , 

where and are respectively the upper and lower sur- 

face stress at the node j of the element i, and a is the 
maximum allowable stress. The nxamber of elements is ne and 
nj is the number of nodal points. 

Displacements : 



u. ./u 
ij max 



1 < 0 



i = l,nj 
j = 1/HC 



where u^^ is the displacement at node i in the coordinate 

direction j and u is the maximiam allowable displacement 

rri3.x 

in coordinate direction j . The number of coordinates is 
nc . 

Height to thickness ratio of frame webs; 



- 5 < 0 



1 - h./t„. < 0 

X Wx — 



i = l,ns 
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are respectively the web height and web 



where h. and t.^. 

1 Wi 

thickness of ith stiffener set. 

Width to thickness ratios of flanges: 

w^/t^^ -5£0 i=l,ns 

1 - w . /t - . < 0 

r f 1 — 

where w^ and t^^ are respectively the width and the thickness 
of flange of the ith stiffener set. 
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IV. DESIGN EXAMPLES 



This chapter will present design examples of unstiffened 
and stiffened pLates. All calculations were carried out on 
the IBM 3033 model 370 computer at the Naval Postgraduate 
School. Material constants used for all examples are Young's 

g 

modulus, E = 3 >: 10 psi, and Poisson's ratio, v = 0.3. Due 
to symmetry only a quarter of the plats was modelled in each 
example . 

A. DESIGN CASE I — THE UNSTIFFENED PLATE 

A quarter of the unstiffened plate simply supported, as 
shown in Fig. 4.1, was modelled using a 4 x4 mesh of the 
plate elements. There are 10 design variables as illustrated 
in Fig. 4.2. The volume of the plate was minimized subject 
to a total of 378 constraints being stress and nodal deflec- 
tion constraints, with a = 20,000 psi and u = 1/2 t. 

max ^ max 

1 . Case lA with A Concentrated Load 

A concentrated load of 1,000 lbs was applied at the 
center of the whole plate. The final optimization results 
show that there are 3 critical deflection constraints of 
the nodes 3, 4 and 10, and that the plate thicknesses of the 
diagonal elements are much thicker than those of off-diagonal 

3 

elements. The final volume is 363.916 in . The results of 
this case are summarized in Table I for comparison with the 
other cases. 
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t : plate thickness 
Q : element 

« ; node 



Figure 4.1 Design Case I — The Unstiffened Plate Model 



36 



t (10) 



t (8 ) 



t (5) 



•:( 1 ) 




Figure 4.2 Design Case I — Design Variables 
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TABLE I 





Summary of Design Case 


I 


'ariable 


Concentrated 


Distributed 




Load 


Load 


t(l) 


1. 302 


0.352 


t (2) 


1.0 31 


0.263 


t (3) 


0,573 


0.589 


t(4) 


1.059 


0.799 


t(5) 


0.120 


0.086 


t{6) 


0.120 


0.085 


t(7) 


0.120 


0.086 


t(8) 


0.120 


0.086 


t (9) 


0.120 


0.086 


t (10) 


0.120 


0.085 


Vol 


363.916 


170.396 


t: 


plate thickness (in) 




Vol ; 


volume (in^) 




2. Case 


IB with Uniform Distributed 


Load 



2 

A vuiiform distributed load, 0.278 Ibs/in was applied 
over the whole plate. The final optimization results show 
that there are 6 critical deflection constraints of the 
nodes 3, 4, 5, 9, 10 and 15, and that the plate thickness 
has a similar trend to case lA where diagonal elements are 
thicker than those of the off-diagonal elements. The final 

3 

volume is 170.396 in . The results are summarized in Table I 
with those of case lA. 

B. DESIGN CASE II — THE STIFFENED PLATE 

A quarter of the stiffened plate simply supported, as 
shown in Fig. 4.3, was modelled 4x4 mesh of the plate ele- 
ments and 12 stiffener elements. 




# : plate/stiffener node 



CASr A 




CASE P 




SECTION A- A AND SECTION B-B 



CASE A 



CASE ? 




SECTIO?: c-c 



Figure 4.3 Design Case II — The Stiffened Plate Model 
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1 . Case IIA — Rectangular Type Frame Stiffener 



There are 17 design variables: a plate thickness, 

and 8 sets of stiffener heights and thicknesses, cis tabu- 
lated in Table II. The volume of the stiffened p].ate was 
minimized subject to a total of 222 constraints; stress, 
stiffener height to thickness ratio, and nodal deflection 
constraints . 

a. Case IIAl with A Concentrated Load 

A concentrated load of 1,000 lbs v/as cipplied at 
the center of the whole plate. There are 14 critical con- 
straints; the height to thickness ratios, and 6 deflections 
of the nodes 4, 5, 10, 29, 30 and 31. The final volume is 

3 

102.017 in with a uniform plate thickness of 0.031 in. The 
results are summarized in Table III for comparing with the 
other cases. 

b. Case IIA2 with Uniform Distributed Load 

A uniform distributed load, 0.278 lb/ in was 
applied over the whole plate. There are 12 critical con- 
straints: the height to thickness ratios, and 4 deflections 

of the plate nodes 7, 8, 14 and 19. The final voliame is 

3 

65.984 in with a uniform plate thickness of 0.059 in. The 
results of this case are summarized in Table III with those 
of the case IIAl. The volumes of the above two cases are 
significantly reduced in comparison to those of the unstif- 
fened plate. 
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TABLE II 



Design 


Case IIA--Design 


Variables and 


Nodal 


Connectivity 


element 

number 


design 

variables 


N1 


node 

N2 


number 

N3 


N4 


pi 


t 


1 


6 


7 


2 


p2 


t 


2 


7 


8 


3 


p3 


t 


3 


8 


9 


4 


p4 


t 


4 


9 


10 


5 


P5 


t 


6 


11 


12 


7 


p6 


t 


7 


12 


13 


8 


P7 


t 


8 


13 


14 


9 


p8 


t 


9 


14 


15 


10 


p9 


t 


11 


16 


17 


12 


plO 


t 


12 


17 


18 


13 


pll 


t 


13 


18 


19 


14 


pl2 


t 


14 


19 


20 


15 


pl3 


t 


16 


21 


22 


17 


pl4 


t 


17 


22 


23 


18 


pl5 


t 


18 


23 


24 


19 


pl6 


t 


19 


24 


25 


20 


si 


t^^(l),h(l) 


30 


35 






s2 


t^-(2) ,h(2) 


35 


36 






s3 


t^.^(3) ,h(3) 


36 


37 






s4 


t^,(4) ,h(4) 


37 


38 






s5 


t^,(5) ,h(5) 


29 


30 






s6 


t^^(6) ,h(6) 


28 


29 






s7 


t^(7) ,h(7) 


27 


28 






s8 


t^(8) ,h(8) 


26 


27 






s9 


t^(5) ,h(5) 


30 


31 






slO 


t^(6),h(6) 


31 


32 






sll 


t^(7),h(7) 


32 


33 






sl2 


t^(8) ,h(8) 


33 


34 
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TABLE III 



Summary of Design Case IIA 



variable 




concentrated 


distribut 






load 


load 


t 




0.051 


0.059 






0.438 


0.209 


t^(2) 




0.318 


0.168 






0.553 


0.267 






0.673 


0.344 






0.097 


0.138 






0.096 


0 .164 






0.095 


0.171 


tw(8) 




0.095 


0.155 


h (1) 




2.189 


1.044 


h (2) 




1.568 


0.832 


h (3) 




2.753 


1.329 


h (4) 




3.343 


1.709 


h (5) 




0.486 


0.687 


h (6) 




0.482 


0.821 


h (7) 




0.478 


0.851 


h (8) 




0.477 


0.776 


Vol 




102.017 


65.984 


t; 


plate 


thickness (in) 




Sj* 


frame 


thickness (in) 




h: 


frame 


height (in) 




Vol: 


volume 


: (in^) 
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2 . Case IIB — T Type Frame Stiffener 



There are 33 design variables: a plate thickness, 8 
sets of the web heights and thicknesses, and 8 sets of the 
flange widths and thicknesses, as tabulated in Table IV. 

The volume of the stiffened plate was minimized subject to 
a total of 230 constraints: stress, web height to thickness 

ratio, flange width to thickness ratio, and nodal deflection 
constraints. 

a. Case IIBl with A Concentrated Load 

A concentrated load of 1,000 lbs was applied at 
the center of the whole plate. There are 12 critical con- 
straints: the web height to thickness ratios, the 3rd and 

8th flange width to thickness ratios, and deflection con- 
straints at the node 5 of the element p5 and node 30 of the 

3 

element si. The final volume is 55.091 in with a uniform 
plate thickness of 0.023 in. The results of this case are 
summarized in Table V for comparison with the other cases. 

b. Case IIB2 with Uniform Distributed Load 

2 

A uniform distributed load, 0.278 Ibs/in was 
applied over the whole plate. There are 16 critical con- 
straints: the web height to thickness ratios except the 

8th one, the 3rd through 8th flange width to thickness ratios 

except the 6th one, and 4 deflections of the plate nodes 7, 

3 

8, 14 and 19. The final voliame is 64.414 in with a uniform 
plate thickness of 0.060 in. This volume is greater than the 
volume of the above concentrated load case, being different 
from the other cases which the latter is greater than the 
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TABLE IV 



Design Case IIB — Design Variables and Nodal Connectivity 



element 


design 




node 


number 




number 


variables 


N1 


N2 


N3 


N4 


pi 


t 


1 


6 


7 


2 


p2 


t 


2 


7 


8 


3 


p3 


t 


3 


8 


9 


4 


p4 


t 


4 


9 


10 


5 


P5 


t 


6 


11 


12 


7 


p6 


t 


7 


12 


13 


8 


P7 


t 


8 


13 


14 


9 


p8 


t 


9 


14 


15 


10 


p9 


t 


11 


16 


17 


12 


plO 


t 


12 


17 


18 


13 


pll 


t 


13 


18 


19 


14 


pl2 


t 


14 


19 


20 


15 


pl3 


t 


16 


21 


22 


17 


pl4 


t 


17 


22 


23 


18 


pl5 


t 


18 


23 


24 


19 


pi 6 


t 


19 


24 


25 


20 


si 


t^(l) ,h(l) ,t^(l) ,w(l) 


30 


35 






s2 


t^(2) ,h(2) ,tj(l) ,w(l) 


35 


36 






s3 


t^,(3) ,h(3) ,tf(3) ,w(3) 


36 


37 






s4 


t^(4) ,h(4) ,tf(4) ,w(4) 


37 


38 






s5 


t^^(5) ,h(5) ,tg(5) ,w(5) 


29 


30 






s6 


t^(6) ,h(6) ,tf (6) ,w(6) 


28 


29 






s7 


t^(7) ,h(7) ,tf (7) ,w(7) 


27 


28 






s8 


CO 

00 

00 

00 

4-> 


26 


27 






s9 


t^(5) ,h(5) ,tf (5) ,w(5) 


30 


31 






slO 


t^(6) ,h(6) ,tg(6) ,w(6) 


31 


32 






sll 


t^(7) ,h(7) ,tf (7) ,w(7) 


32 


33 






si 2 


t^(8) ,h(8) ,t^(8) ,w(8) 


33 


34 
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TABLE V 



Summary of Design Case IIB 



variable 


concentrated 


distribut 




load 


load 


t 


0.023 


0.060 


t^(l) 


0.133 


0.14 9 




0.077 


0.126 


t^(3) 


0.137 


0.197 


t^(4) 


0.162 


0.242 




0.381 


0.103 


tw(6) 


0.355 


0.114 




0.316 


0.115 


t^(8) 


0.274 


0.122 


h (1) 


0.659 


0.736 


h (2) 


0.386 


0.631 


h (3) 


0.680 


0.976 


h (4) 


0.805 


1.201 


h (5) 


1.902 


0.509 


h (6) 


1.767 


0.566 


h (7) 


1.570 


0.569 


h (8) 


1.368 


0.602 


1 — 1 
4-1 


0.021 


0.262 


tf (2) 


0.028 


0.183 


tfO) 


0.011 


0.308 


tj(4) 


0.026 


0.460 


tf (5) 


0.876 


0.209 


tg(6) 


0.799 


0.231 


tf (7) 


0.6 71 


0.257 


tj(8) 


0.460 


0.201 
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TABLE V (Cont’d) 



variable 


concentrated 


di.stributi 






load 


load 


w (1) 




0.083 


0.278 


w (2) 




0.054 


0.187 


w (3) 




0.057 


0.307 


w (4) 




0.074 


0.460 


w (5) 




0.891 


0.210 


w (6) 




0.814 


0.235 


w (7) 




0.685 


0.258 


w (8) 




0.464 


0.201 


Vol 




55.091 


64.414 




t : 


place thickness (in) 






W 


web thickness (in) 






h: 


web height (in) 






tf : 


flange thickness (in) 






w: 


flange width (in) 





Vol: volume (in-’) 
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former. This suggests that the stiffener type and arrange- 
ment used here is particularly efficient for the concentrated 
load case. The results of this case are sxommarized in Table 
V with those of the case IIBl. The volumes of the above 
cases are significantly reduced with comparison to those of 
unstiffened plates. The volume of T type stiffener case 
is slightly smaller than that of rectangular case in the 
distributed load. On the other hand, the volume of T type 
stiffener case is much smaller than the rectangular case 
for the concentrated load. 

C. DESIGN CASE III — THE HATCH COVER 

A quarter of a hatch cover was modelled 4x4 mesh of the 
plate elements and 24 stiffener elements, as shown in Fig. 

4.4. 

1 . Case IIIA-- Rectangular Type Frame Stiffener 

There are 9 design variables: a plate thickness, 

and 4 sets of stiffener heights and thicknesses, as tabulated 
in Table VI. The volume of the hatch cover was minimized 
subject to a total of 274 constraints; stress, height to 
thickness ratio, and nodal deflection constraints, 
a. Case IIIAl with A Concentrated Load 

A concentrated load of 1,000 lbs was applied at 
the center of the whole plate. There are 4 critical con- 
straints; the 2nd and 4th height to thickness ratios, and 
the upper and lower stresses of the node 5 of the element 
p4 . The dimensions of the 2nd stiffener set are negligibly 
small with comparison to those of the other stiffener sets. 
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Figure 4.4 Design Case III — The Hatch Cover Model 
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TABLE VI 



Design Case IIIA — Design Variables and Nodal Connectivity 



element 


design 




node 


number 




number 


variables 


N1 


N2 


N3 


N4 


pi 


t 


1 


6 


7 


2 


p2 


t 


2 


7 


8 


3 


p3 


t 


3 


8 


9 


4 


p4 


t 


4 


9 


10 


5 


p5 


t 


6 


11 


12 


7 


p6 


t 


7 


12 


13 


8 


P7 


t 


8 


13 


14 


9 


p8 


t 


9 


14 


15 


10 


p9 


t 


11 


16 


17 


12 


plO 


t 


12 


17 


18 


13 


pll 


t 


13 


18 


19 


14 


pi 2 


t 


14 


19 


20 


15 


pi 3 


t 


16 


21 


22 


17 


pi 4 


t 


17 


22 


23 


18 


pi 5 


t 


18 


23 


24 


19 


pi 6 


t 


19 


24 


25 


20 


si 


l/2t^(3) ,h(3) 


26 


27 






s2 


l/2t^(3) ,h(3) 


27 


28 






s3 


t^(l) ,h(l) 


28 


29 






s4 


t^(l) ,h(l) 


29 


30 






s5 


tw(4),h(4) 


26 


31 






s6 


t^(2) ,h(2) 


28 


32 






s7 


t^(l) ,h(l) 


30 


33 






s8 


t^(4) ,h(4) 


31 


34 






s9 


t^(2) ,h(2) 


32 


36 






slO 


t^(l) ,h(l) 


33 


38 






sll 


t^,(3) ,h(3) 


34 


35 






sl2 


t^(3) ,h(3) 


35 


36 
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TABLE VI (Cont'd) 



element 



design. 



nxamber 


variables 


N1 


N2 


si 3 


t^^,(2),h(2) 


36 


37 


sl4 


t./2),h(2) 


37 


38 


si 5 


t^^(4),h(4) 


34 


39 


sl6 


t.^(3) ,h(3) 


36 


40 


si 7 


l/2t^,(3) ,h(3) 


38 


41 


sl8 


t^(4) ,h(4) 


39 


42 


sl9 


t^{3) ,h(3) 


40 


44 


s20 


l/2t^,(3) ,h(3) 


41 


46 


s21 


t^^.{4) ,h(4) 


42 


43 


s22 


t^j{4) ,h(4) 


43 


44 


s2 3 


t^,(4) ,h(4) 


44 


45 


s24 


t^^,(4),h(4) 


45 


46 



node number 
N3 



N4 
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This suggests that the 1st stiffener set is stiff enough to 
compensate the elements attached to the 2nd stiffener set. 

3 

The final volume is 323.639 in with a uniform plate thick- 
ness of 0.359 in. The results of this case are summarized 
in Table VII. 



TABLE VII 

Summary of Design Case IIIA 



variable 


concentrated 

load 


distributed 

load 


t 


0.359 


0.059 


ft 


0.811 


0.863 




0.001 


0.002 




0.084 


0.518 




0.015 


0.009 


h (1) 


1.306 


1.750 


h (2) 


0.006 


0.011 


h (3) 


0.162 


1.484 


h (4) 


0.073 


0.019 


Vol 


323.639 


87.314 



t: plate thickness (in) 

t^: frame thickness (in) 

h: frame height (in) 

Vol; volume (in^) 



b. Case IIIA2 with Uniform Distributed Load 

2 

A uniform distributed load, 0.278 Ibs/in was 
applied over the whole plate. There are 6 critical constraints 
the 2nd height to thickness ratio, and 5 deflections of the 
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plate nodes 8, 9 and 14 and of stiffener nodes 32 and 37. 

With the same suggestion as the concentrated load case the 
dimensions of the 2nd stiffener set are negligibly small. 

3 

The final volume is 87.314 in with a uniform plate thickness 
of 0.059 in. The design results of the above both cases 
are summarized in Table VII with those of case IIIAl. The 
volumes of the above two cases are significantly reduced 
with comparison to those of unstiffened plates. 

2 . Case IIIB — T Type Frame Stiffener 

There are 17 design variables: a plate thickness, 4 

sets of the web heights and thicknesses, and 4 sets of the 
flange widths and thicknesses, as tabulated in Table VIII. 

The volume of the hatch cover was minimized subject to a t^tal 
of 278 constraints: stress, web height to thickness ratio, 

flange width to thickness ratio, and nodal deflection 
constraints . 

a. Case IIIBl with A Concentrated Load 

A concentrated load of 1,000 lbs was appi;Lsd at 
the center of the whole plate. There are 3 critical con- 
straints: the 3rd flange width to thickness ratios, and 2 

deflections of the plate node 5 and stiffener node 30. 
Negligibly small are the whole dimensions of the 2nd stiffener 
set and the flange dimensions of the 3rd stiffener set with 

comparison to the other dimensions. The final volume is 
3 

95.883 in with a uniform plate thickness, 0.083 in. 
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TABLE VIII 



Design 


Case IIIB--Design Variables 


and 1 


[^odal 


Connectivity 


element 


design 




node 


number 




number 


var icibles 


N1 


N2 


N3 


N4 


pi 


t 


1 


6 


7 


2 


p2 


t 


2 


7 


8 


3 


P3 


t 


3 


8 


9 


4 


p4 


t 


4 


9 


10 


5 


p5 


u 


6 


11 


12 


n 


p6 


t 


7 


12 


13 


8 


pv 


t 


8 


13 


14 


9 


p8 


t 


9 


14 


15 


10 


p9 


t 


11 


16 


17 


12 


plO 


t 


12 


17 


18 


13 


pll 


t 


13 


18 


19 


14 


pi 2 


t 


14 


19 


20 


15 


pl3 


t 


16 


21 


22 


17 


pl4 


t 


17 


22 


23 


18 


pl5 




18 


23 


24 


19 


pi 6 


t 


19 


24 


25 


20 


si 


l/2t^(3) ,h(3) ,t^{3) ,l/2w(3) 


26 


27 






s2 


l/2t^(3) ,h(3) ,tj{3) ,l/2w(3) 


27 


28 






s3 


t^^(l) ,h(l) ,tg(l) ,w(l) 


28 


29 






s4 


t^,(l) ,h(l) ,t^(l) ,w(l) 


29 


30 






s5 


t^(4) ,h(4) ,tf (4) ,w(4) 


26 


31 






s6 


t^(2) ,h(2) , tf (2) ,w(2) 


28 


32 






s7 


t.^(l) ,h(l) ,tf (1) ,w(l) 


30 


33 






s8 


t^(4) ,h(4) ,tf (4) ,w(4) 


31 


34 






s9 


t^(2) ,h(2) (2) ,w(2) 


32 


36 






slO 


t^(l) ,h(l) ,tf (1) ,W{1) 


33 


38 






sll 


t^(3) ,h(3) ,tf (3) ,w{3) 


34 


35 






sl2 


t^(3) ,h(3) ,tf (3) ,w{3) 


35 


36 







53 



TABLE VIII (Cont'd) 



element 


design 




node 


number 


variables 


N1 


N2 


sl3 


t^(2) ,h(2) ,t^(2) ,w(2) 


36 


37 


sl4 


t^(2) ,h(2) ,t^(2) ,w(2) 


37 


38 


sl5 


t^(4) ,h(4) ,tf(4) ,w(4) 


34 


39 


sl6 


t^(3) ,h(3) (3) ,w(3) 


36 


40 


sl7 


l/2t^(3) ,h(3) ,tj(3) ,l/2w(3) 


38 


41 


sl8 


t^(4) ,h(4) ,tf (4) ,w(4) 


39 


42 


sl9 


t^(3) ,h(3) (3) ,w(3) 


40 


44 


s20 


l/2t^(3) ,h(3) (3) ,l/2w{3) 


41 


46 


s21 


t^(4) ,h(4) ,tj(4) ,w(4) 


42 


43 


s22 


t^^(4) ,h(4) ,tf (4) ,w(4) 


43 


44 


s23 


t^(4) ,h(4) ,tf (4) ,w(4) 


44 


45 


s24 


t^^(4) ,h(4) (4) ,w(4) 


45 


46 



number 
N3 N4 
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b. Case IIIB2 with Uniform Distributed Load 

2 

A uniform distributed load, 0.278 Ibs/in was 
applied over the whole plate. There are 9 critical con- 
straints; the 3rd web height to thickness ratio, the 3rd 
flange width to thickness ratio, 4 deflections of the plate 
nodes 8, 9, 13 and 14, and 3 deflections of the stiffener 
nodes 32, 36 and 37. Negligibly small are the v/hole dimen- 
sions of the 2nd stiffener set, and the flange dimensions o 
the 3rd and 4th stiffener sets with comparison to the other 
dimensions. This suggests that the 1st stiffener set is 
stiff enough to compensate the elements attached to the 2nd 
stiffener set. Additionally it can be suggested that the 
3rd and 4th stiffener sets are stiff enough to support the 
given load with only the rectangular type ones. The final 
volume is 76.436 in with a uniform plate thickness, 0,062 
in. The results of the above both cases are summarized in 
Table IX. The volumes of the above two cases are signifi- 
cantly reduced with comparison to those of the unstiffened 
plates and the volxme of T type stiffener case is much 
smaller than that of the rectangular type case. 
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TABLE IX 



Summary of Design Case IIIB 



variable 


concentrated 




load 


t 


0.083 


f — 1 
> 


0.284 




0.003 




0.32 3 


t^(4) 


0.020 


h (1) 


1.201 


h (2) 


0.010 


h (3) 


1.4 33 


h (4) 


0.060 


tg(l) 


1.001 


tf (2) 


0.003 


t^{3) 


0.000 


tf (4) 


0.010 


w (1) 


1.110 


w (2) 


0.010 


w (3) 


0.002 


w (4) 


0.0 30 


Vol 


95.383 


t 


plate thickness (in) 




web thickness (in) 


h 


web height (in) 


t .. 
r 


flange thickness (in 


w 


flange width (in) 


Vol 


volume (in^) 



distributed 

load 

0.062 

0.284 

0.003 

0.318 

0.020 

1.200 

0.010 

1.430 

0.060 

1.000 

0.003 

0.001 

0.010 

1.110 

0.010 

0.003 

0.001 

76.436 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

A finite element analysis program was developed and 
coupled to an optimizer for stiffened plate design. Stif- 
fened plates were designed for minimum volume with a concen- 
trated load. 

Five important conclusions are made from considering 
several design examples. 

First, a stiffened plate is much more efficient than 
the unstiffened plate for the optimum volume design. 

Second, it is required to attach the larger stiffeners 
to the diagonal elements of the plate as compared to the 
off-diagonal elements. 

Third, the T type stiffeners provide a lower volume 
design than the rectangular type stiffeners. 

Fourth, the stiffened plate with the diagonally attached 
stiffeners are much more efficient than those with the 
rectangularly attached stiffeners. 

Fifth, the stiffened plate attached diagonally with the 
T type stiffeners is particularly efficient for the concen- 
trated load case. 

These conclusions are based on the design of simply 
supported plates and the conclusions could be quite differ- 
ent if other boundary conditions are considered. 
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B. RECOMMENDATIONS 



The following recommendations may be useful for future 
work . 

1. The analysis program should be modified to a 
generalized program with a variety of elements. 

2 . The analysis program should be extended to 
consider the dynamic problem. 

3. The analysis program should be extended to 
consider the multiple loading conditions. 

4 . Routines should be added to calculate gradients 
analytically. 

5. Buckling constraints of the plate and stiffener 
elements should be added in the optimization 
process . 

6 . Frequency constraints should be added in the 
optimization process. 

7. The analysis program should be modified to reduce 
the computing time during the optimization process. 
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